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Abstract 

We propose a new method for the nonperturbative solution of quantum held theories and illus- 
trate its use in the context of a light-front analog to the Greenberg-Schweber model. The method 
is based on light-front quantization and uses the exponential-operator technique of the many-body 
coupled-cluster method. The formulation produces an effective Hamiltonian eigenvalue problem in 
the valence Fock sector of the system of interest, combined with nonlinear integral equations to be 
solved for the functions that define the effective Hamiltonian. The method avoids the Fock-space 
truncations usually used in nonperturbative light-front Hamiltonian methods and, therefore, does 
not suffer from the spectator dependence, Fock-sector dependence, and uncanceled divergences 
caused by such truncations. 
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I. INTRODUCTION 

The central problem of a quantum field theory is to compute its mass spectrum and the 
corresponding eigenstates. All physical quantities can be computed from these. If the theory 
is quantized in terms of light-front coordinates [H, this spectral problem can be written as 
a Hamiltonian eigenvalue problem [2|, V^liplP)) = P^\ip{P))i where = — is 
the light-front energy operator, V_ = (P"*" = P° + V^,V± = {V^,V^)) is the light-front 
momentum operator, and P'^ are the corresponding eigenvalues. For an eigenstate of mass 
M, the mass-shell condition = yields P~ = {M'^ + Pl)/P+. Thus, eigenvalues of V~ 
determine the mass spectrum. 

The standard light-front Hamiltonian approach is to expand {ipiP)) in a set of Fock states, 
eigenstates of P with definite numbers of constituents. The coefficients in the expansion are 
the light-front momentum-space wave functions. This takes advantage of two important 
aspects of light-front coordinates the relative-momentum coordinates separate from the 
external momentum, so that the wave functions depend only on the relative momenta, and 

the positivity of P~^ = \J (P)^ + + P^ excludes vacuum contributions to the expansion, 
so that the wave functions represent the properties of the eigenstate only. 

Given the Fock-state expansion, the eigenvalue problem becomes an infinite set of cou- 
pled integral equations for the wave functions. The expansion and the coupled system are 
truncated to yield a finite problem, which is then solved, usually by numerical techniques 

In more than two dimensions, some form of regularization is required to properly define 
the integrals of the coupled system. The cancellations that must take place in the regu- 
larization scheme are disrupted by the truncation, resulting in uncanceled divergences. A 
re-parameterization of the theory, such as sector-dependent parameterization |#^, can be 
arranged to appear to absorb these divergences, but not simultaneously for all physical quan- 
tities [8!]. The truncation also causes self-energy contributions and vertex functions to be 
dependent on the momenta of Fock-state constituents that are only spectators to the process 
in question. This spectator and Fock-state dependence results in great complications for the 
analysis and solution of the theory. 

In particular, the Ward identity of gauge theories is destroyed by truncation. For photon 
emission in QED, a one-photon truncation keeps only the self-energy correction to the elec- 
tron leg on the side opposite the photon emission; the self-energy correction on the other 
leg and the vertex loop correction are eliminated js], lof. The relevant diagrams are shown in 
Fig. [T| only the first survives the truncation. Thus, the Ward identity connecting vertex and 



FIG. 1. Graphs contributing to the Ward identity in QED. Only the first contributes in a one- 
photon truncation of the Fock space. 

wave-function renormalization is broken for interactions internal to a bound-state problem. 
This is what drives the renormalization of the charge in a sector-dependent parameterization 
of the theory [s], , but this is clearly unphysical and has nothing to do with ordinary charge 
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renormalizat ion 

The analog of these difficulties with truncation can be induced in Feynman perturbation 
theory by separating covariant diagrams into time-ordered diagrams and discarding those 
time orderings that include intermediate states with more particles than some finite limit. 
This destroys covariance, disrupts regularization, and induces spectator dependence for sub- 
diagrams. In the nonperturbative case, this happens not just to some finite order in the 
coupling but to all orders. 



II. LIGHT-FRONT COUPLED-CLUSTER METHOD 

To avoid truncation, we introduce the exponential-operator technique of the many-body 
coupled-cluster (CC) method H, 13|. An eigenstate \ip{P_)) is written as y/Ze'^\(t>{P) , where 



|0(P)) is limited to one or a few Fock sectors with the lowest number (s) of constituents, the 
valence sector (s). The operator T is a sum of operators that increase particle number but 
conserve momentum P, the angular momentum componenlH J^, and all relevant quantum 
numbers, such as charge and baryon number. The factor y/Z is a normalization factor, such 
that = = ^{E^ - P)- We then construct an effective eigenvalue 

problem in the valence sector, PuP~|0(P)) = — p+ ^ where V~ = e V~e and P„ is 

a projection onto the valence sector. Equations for the functions that determine T are found 
by the orthogonal projection (1 — Pt,)P^|0(P)) = 0. Up to this point, no approximation has 
been made, and the problem remains infinite, because there are infinitely many contributions 
to T. 

To have a finite set of equations, we truncate T to a few operators and truncate the 
projection 1 — P^ in a consistent way, such that just enough equations are produced to be 
able to solve for the functions in the truncated T operator. For example, if T can create 
one additional particle above the valence state, 1 — P„ projects onto only this additional 
Fock sector. After truncation, we have a finite set of nonlinear equations for the functions 
in T, coupled to the valence-sector wave functions, and a valence-sector eigenvalue problem 
where the effective Hamiltonian depends on the functions in T. The former are essentially 
auxiliary equations that help define the latter. 

What is not truncated is the exponentiation of T, and thus the full Fock space can be 
retained, though the wave functions for the higher Fock sectors are clearly only approximate. 
The effective Hamiltonian is computed from its Baker-Hausdorff expansion V~ = V~ + 
[P^, T] + |[[P^, T], T] + • ■ ■ . Only a finite number of terms contributes, because each factor 
of T increases the number of particles created, eventually exceeding the truncation of (1— P^,). 

Although this light-front coupled-cluster (LFCC) method uses the mathematics of the 
traditional CC method [iMl, it is quite different conceptually. In fact, the name coupled 



cluster does not really apply, but we use it to acknowledge the origin of the LFCC method. 
The CC method is applied to a single Fock sector, with a large number of constituents. The 
T operator builds correlated excitations onto a Hartree-Fock-type ground state. Within 
products of T there are no contractions, because every term in T annihilates one or more 
of the single-particle states in the ground state and creates one or more excited states. In 
the LFCC method, the valence sector has a small number of constituents, and the method 

^ For external photon emission, the truncation does not apply and the Ward identity is preserved. Without 



vacuum polarization, the plus component of the dressed-electron current is not renormalized 10|, 111] . 



^ The other two components of angular momentum are not kinematic ^ . The eigenstates of the Hamiltonian 
will in general be linear combinations of eigenstates of ,P . Determination of the eigenstates is a separate 
dynamical problem. 



of solution of the eigenvalue problem here is left unspecified. The terms of the T operator 
do include annihilation, because the positive light-front momentum P+ cannot be conserved 
unless one or more particles are annihilated to provide momentum for those that are created. 
As a consequence, powers of T include contractions, but these are needed in order that T 
to some power not annihilate the entire valence state, which would effectively truncate the 
exponentiation of T. 

In addition to the fundamental mass eigenvalue problem, the LFCC method must also 
contend with the evaluation of matrix elements of operators, in order to be able to extract 
physical quantities from the LFCC eigenstates. This is nontrivial, because a direct calcu- 
lation of the normalizing factor \fZ is impractical, due to the infinite set of terms in the 
sum over Fock states within (0(£)|e^^e^|0(P)). This same issue arises in the traditional CC 



method [13|], and there a technique exists for expectation values which can be adapted for 
the LFCC method and extended to include off-diagonal matrix elements. Some care must 
be taken, however, in that the LFCC method uses momentum eigenstates with Dirac-delta 
normalization, unlike the unit normalization of the standard CC states. The normalization 
factor \fZ is introduced to avoid division by (V^(P')lV^(i^)) = ^{El ~ E) in the computation 
of expectation values. 

For an operator O we write the expectation value (O) in the state \/Ze^\(f){P_)) as (O) = 
Z((/)(P)|e^^Oe^|0(P)) and define O = e-^Oe^ and 



(^(P)| = Z{<j,{P)\e^'e^ = v^(^(P)|e^, (2.1) 
so that (O) = {i^{P.)\0\(f){P_)). By construction, we have 



(v^(£')i0(£)) = mEmp)) = 5{E - p) (2.2) 

and 

(^(P)l^- = v^(^(P)|e^e-^p-e^ = _^(^(P)|. (2.3) 

Thus, {ip{P) \ is a left eigenvector of the (necessarily) non-Hermitian , with the same mass 
eigenvalue, normalized such that the projection onto the valence state is a simple momentum- 
conserving delta function. Therefore, an expectation value is computed by constructing the 
effective operator O from a Baker-Hausdorff expansion, solving the left-hand eigenvalue 
problem, and evaluating the inner product ('?/'(P)|O|0(P)). As for V~ , only a finite number 
of terms in the Baker-Hausdorff expansion of O will contribute. The extension to off-diagonal 
matrix elements is straightforward. 

The left-hand eigenvalue problem must be truncated to an extent consistent with the 
truncation of T, such that {'ip{P)\ is limited to the Fock sectors of the valence state plus 
those created by application of T. To understand the truncation, consider the following. 
Define an operator L = {l — P.u)Ze^''e^P^^ with P^ the projection onto the valence eigenstate 
|0(£))- Because of the projection operators, is simply 1 + L. The left-hand eigenstate 
can then be written as (V^(P)| = (0(P)|e^^ + Z(0(P)|e^^e^(P^ - P^y . We see, then, that 
L plays the role of T, and therefore should be truncated in the same way. The truncated 
left-hand eigenvalue problem creates a finite set of linear equations for the wave functions 
of {^{P)\. 
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III. MODEL APPLICATION 



To illustrate the method, we apply it to a simple model where an analytic solution is 



known [15|. The model is a light-front analog of the Greenberg-Schweber model [16[ for a 



static fermionic source that emits and absorbs bosons without changing its spin. In updated 



ai{q) (3.1) 



notation, the Hamiltonian given in [15| can be written as 

where a\ creates a "physical" boson of mass yUo, o\ creates a Pauli-Villars (PV) boson of 
mass /ii, and h\ creates the fermion with mass M and spin s. The parameter 7 can take 
any positive value; as shown in [isj, it controls the longitudinal endpoint behavior of the 
wave functions. The PV boson provides the necessary ultraviolet regularization, to define 
the self-energy M'q. To accomplish the regularization, the PV boson is assigned a negative 
normJl The (anti) commutation relations are 

{K{v)A'ii)} = Sss'Sip-^), h(g),4(^')] = {-iySu'S{q-q;). (3.2) 
A graphical representation of the Hamiltonian is given in Fig. |2J The model is not fully 



« — X — * — - ^ ^. — — .-^ — . — — » r « — 

FIG. 2. Graphical representation of the model Hamiltonian operator V~ defined in Eq. (|3.ip of 
the text. Each graph represents an operator that annihilates one or more particles on the right 
and creates one or more to take their place. The crosses refer to hght-front kinetic-energy terms. 

covariant, which hides some of the power of the LFCC method, but is sufficient to show how 
the method can be applied. 

We truncate the T operator to include only boson emission from the fermion, as repre- 
sented in Fig. 131 

T^^f^MMmilMl^l), (3.3) 

Is 

with tis the operator functions to be determined. The effective Hamiltonian V~ is con- 



< ^ 



FIG. 3. Graphical representation of the truncated T operator. 



In [15| , the PV cancellations were arranged by use of an imaginary coupling rather than a negative norm. 
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structed from its Baker-Hausdorff expansion. A graphical representation of the first two 
commutators is given in Fig. |H The expression for V~ is 

,2 _L ^2 

(3.4) 

+ M'qP+ 



P+ 
dq 



dpdq 



wq 

+ \ 7 



-yi ^ 

V/I67r3g+ ^ \p+ + g+ 
c/pc/g ^ a\ (q) hi (p) bs {p + q) 



blip + q)bs{p)ai{q) 



Is 



X 



P+ V/I67r3g+ VP"^ + q 



2P+ 



+ 



/if + ql M'^q+ 
q+ P+ 



tis{q,p) 



pi 



q'^)tis{q, p - (/)ti's{q^, p + q-<f)} 



-29 {p-^ + g+ - g'+ 
dpdqdq' 



p + q — q 



tu{q,p)tv,{q\p + q - q') 



P+ 



6'(p^ — g 



P+ _|_ gH 

' - ^'^) 5Z + q- i)bs{p)ai{q) 

P+ _ g' 



ll's 

\ 7 



g+ _ q'+ 



ti's{q',p- q') 



+ , + , ti's{q',P + q-q') 
^p+ -\- q+ J - - - - _ 

where we hst only terms that connect the lowest Fock sectors. Notice that the self-energy 
contribution Mq is the same in all Fock sectors and that Fig. Hl^b) contains all three of the 
diagrams analogous to those for the Ward identity in QED, as discussed in the Introduction, 
with no truncation in particle number. 

The valence state is the bare-fermion state \4>"{P_)) = bl^{P_)\0). The projection 1 — 
Py is truncated to the one-fermion/one-boson sector. The full eigenstate is \ip'^{P_)) = 
a/Zc^ 1 (P) ) , where we have generalized the basic construction to have one T operator for 
both spins a = ± and will solve for both states simultaneously. This allows contractions in 
powers of T to include sums over both spins. The truncated left-hand eigenvector is 



Is 



(3.5) 



where is the left-hand one-fermion/one-boson wave function. Due to the lack of covariance 
in the model, these states are all limited to having a fixed total transverse momentum P±, 
which we take to be zero. 
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(a) 



■ t ■ y 



(b) 



FIG. 4. Graphical representation of the operators (a) [7^~,T] and (b) [['P~,T],T]. The crosses 
indicate hght-front kinetic-energy contributions. The self-energy loops in the fourth diagram of 
(a) and the first and second diagrams of (b) make identical contributions, with no Fock-sector or 
spectator dependence. 



The eigenvalue problem in the valence sector P^V \(t)'^{P)) = prl^'^d^)) becomes 
M2 + M^P+ , 



+ 



dq 



9 

P+ 



-.9{P+ - 



P+ -q 
P+ 



+ \ 7 



5^(-i)V(g,£-g)|0^(a) 



which reduces to a determination of the self-energy 



'I67r3g+ ^ ' ' V ^ 

In the one-fermion/one-boson sector, we have (1 — Pv)'P~\'P'^{P)) = or 
dq9{P+ ~ q+)Y,A{q)bUP- qm 

g 1 ( r' -q 



X 



v/l67r3g+ 
dq' 



9 



2P^ 



P+-g+-g'+\^ 



X {t/±(g,£-g)t/'±(g',P-g-g') 

+ ^(P+ - (g, P - g - g')ip± (^', £ " } 

-2^(P+ - g'+ 



P+ -g'+\^ 



P^ 



ti±{q,P-q)ti'±{q',P-q') 



(3.6) 



(3.7) 



(3.8) 
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Thus, the contents of the outer curly brackets must sum to zero. This will occur if the 
function tis is 



-9 



p 

P+ _|_ gH 



(3.9) 



The fact that the self-energy Mq is the same in the valence sector and the one-fermion/one- 
boson sector plays a critical role; the expression (13. 7p obtained in the valence sector is 
exactly what is needed to obtain the necessary cancellations in (13. 8p . The self-energy can 
be computed from Eq. (13. 7p as 



ln(/ii//io) 



167r3p+ 7 + 1/2 



(3.10) 



which agrees with the result in fld\. In fact, with tis as given above, the exponential operator 
generates the exact solution given in 15 . 



The solution for tis provides the input to the left-hand eigenvalue problem, {ip'^iE^lV^ = 
^(V'^(£)|- The effective Hamiltonian V~ simplifies considerably; the square bracket in the 
b'^b term becomes just M'^ / and the entire a^Wb term, which corresponds to the curly 
brackets in (13.80 . is zero. The remaining terms in V~ yield the following integral equation 
for the left-hand wave function: 
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1 



_9_ 



2 



P+ 



Ss 



(3.11) 



dq' 



^16n^q+ 



e{p' 



q' -q 



q^ 



(P 



167r3 



— q 



I' 



fif, + q'l 



(P^ 
P+ 



g/+)7(P+ - q 



+ 17 



'+ \ 7' 



0. 



p+ J \ p 

Following the pattern of the inhomogeneous term, we can seek a solution of the form 

-9 fP 



lis{q,P) = Sa 



16^rV V P 



(3.12) 



Substitution yields a one- dimensional integral equation for /(y) 



l{y) = 1 + 



9' /i? 



167r2 



dy'{l-y'y'y'[{l-yyi{y'{l-y))-l{y')]. 



(3.13) 



This equation can be solved iteratively, to generate an expansion in powers of g"^, or numer- 
ically. A Gauss- Jacobi quadrature will convert the integral equation into a linear system 
for the values of / at the chosen quadrature points. The solution then provides the rest of 
the information needed for the computation of matrix elements. 

To consider a particular matrix element as an example, we compute the Dirac form 
factor for the dressed fermion from a matrix element of the current = ip'^^ip. The 
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current couples to a photon of momentum q. With our normalization, the matrix element 
is generally [l3| 

(V^-(P + g)|16vrV+(0)|^±(P)) = 25^^F,{q') ± l^5^^F,{q'), (3.14) 

with Pi and P2 the Dirac and Pauli form factors. In the present model, the fermion cannot 
flip its spin; therefore, P2 is zero, and we investigate only Pi. 

In the LFCC method, the form factor is given by the matrix element 

Pi(g2) = M'{i^^{P + q)\jW)\^^m). (3.15) 



with J+(0) = J^(0) + [J+(0),T] + For this model, there are no contributions from 
fermion-antifermion pairs, so that 

and only the flrst two terms of the Baker-Hausdorff expansion contribute to the matrix 
element. The second term is 

dp' f dp 



[J+(0),T] = 2^/ 7=! / 7=f / (3-17) 
The flrst term contributes l/Sn^ to the matrix element; the second contributes 



(V^±(P + g)|[J+(O),r]|0±(P)) = -ij5^(-l)' f dq;e{P+ + q+ - q'^) 

^ I 

Because the model limits calculations to a flxed total transverse momentum, we calculate 
the matrix element in a frame where q± is zero and g"*" is notEl With a = q^ jP^ and 
P' = P + q, we have 

g2 = (p' _ p f = 2M^ - P'+P- - P'-P+ 

= 2M^ - M^il + a) = . 3.19 

On substitution of the solutions for the wave functions and evaluation of the transverse 
integral, the form factor can be written as a function of a 



\yl{y)y{l-yy[l-{l + a)y\ 



- / dyl{y)y{l-yf^ 



(3.20) 

The PV dependence is easily removed in the limit of an inflnite PV mass (/ii — t- 00). If I 
is computed in quadrature, the integrals remaining in Pi can be computed from the same 
quadrature rule for any chosen value of a. If / is instead constructed as an expansion in 
(7^, Pi can also be constructed as an expansion. In any case, in the limit of — ?■ 0, we 
have from (13.191) that a = and, because the two integrals in (13.201) then cancel, Pi(0) = 1, 
consistent with the unit charge in the current J"*" = ip'j'^ip. 

^ In [1^ the matrix element was computed in a frame where = 0. This could be done because the wave 
functions of the exact solution were taken to be boost invariant. Here, although the right-hand eigenvalue 
problem has accidentally provided the exact solution, we continue with the LFCC approximation in the 
calculation of the matrix element, to provide a mgre complete illustration of the method. 



IV. SUMMARY 



We have proposed a new Hamiltonian method for the nonperturbative solution of quan- 
tum field theories that avoids Fock-space truncations. The full eigenstate is constructed 
from the action of an exponentiated operator T on a valence state This yields a va- 
lence eigenvalue problem = ^ with = e~^V~e^, and a set of auxiliary 
equations for the functions in T. Expectation values are computed as (O) = {ip\e~'^Oe'^\(f>), 
with use of the left-hand eigenstate The method then generates approximations by 
truncation of T rather than of Fock space. 

The application to the simple model in Sec. IIIII shows that the construction of gen- 
erates self-energy contributions that are Fock-sector and spectator independent. Thus, the 
uncanceled divergences that can arise from Fock-space truncations do not occur; the self- 
energy is the same in every sector. The application also shows that a simple approximation 
for the T operator can provide a very good approximation to the eigenstate; in this special 
case, the eigenstate is exact. The calculation of a matrix element is demonstrated in the 
calculation of a Dirac form factor. 

The LFCC method is not limited to any particular theory or model, nor to Pauli-Villars 
regularization. It should be applicable to any regularized field theory. Work on an applica- 
tion to QED is in progress, with some preliminary discussion given in [l8|]. For theories with 
symmetry breaking and vacuum structure, modes of zero longitudinal momentum [l9[ play 
some role and would require extension of the method to include them; in particular, the con- 
tributions to the T operator would not require annihilation, and the exponentiation would 
produce generalized coherent states. For discrete light-cone quantization (DLCQ) j^, |20| . 
where longitudinal momentum fractions are restricted to integer multiples of a fundamen- 
tal amount 1/K, truncation to K particles is automatic; however, the method could still 
be applied as a way of reducing the effective dimension of the underlying matrix eigen- 
value problem and allowing higher resolution. Even the supersymmetric form of DLCQ 
(SDLCQ) [111 should be amenable; instead of constructing V~ from a discretized super- 
charge Q~ via = {Q^ ,Q^}/2\/2, to retain the supersymmetric spectrum, the effective 
Hamiltonian would be constructed from effective supercharges Q~ = e~'^Q~e^ . Thus, 
there is a wide range of applications to consider. 
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